Doping Evolution of Oxygen K-edge X-ray Absorption Spectra in Cuprate 

Superconductors 
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We study oxygen A"-edge x-ray absorption spectroscopy (XAS) and investigate the validity of the 
Zhang-Rice singlet (ZRS) picture in overdoped cuprate superconductors. Using large-scale exact 
diagonalization of the three-orbital Hubbard model, we observe the effect of strong correlations 
manifesting in a dynamical spectral weight transfer from the upper Hubbard band to the ZRS band. 
The quantitative agreement between theory and experiment highlights an additional spectral weight 
reshuffling due to core-hole interaction. Our results confirm the important correlated nature of the 
cuprates and elucidate the changing orbital character of the low-energy quasi-particles, but also 
demonstrate the continued relevance of the ZRS even in the overdoped region. 
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Introduction- Unraveling the nature of low-energy 
quasi-particles in cuprate superconductors is crucial 
to understanding their unconventional superconducting 
mechanism. Despite two decades of studies, however, 
the minimal model for describing the quasi-particle band 
which emerges upon doping remains unclear. Previous 
theoretical work has focused on the single-band Hub- 
bard model[l] or the t — J model obtained by project- 
ing out doubly-occupied charge configurations. [2, 3] In 
these down-folded Hamiltonians, the fundamental quasi- 
particle has been assigned to the so-called Zhang-Rice 
singlet (ZRS):[3, 4] a locally bound d 9 copper 3d x 2_ y 2 
hole hybridized with a doped ligand hole (L) distributed 
on the planar oxygen 2p XyV orbitals [Fig. 1(a)]. 

The relevance of the ZRS in cuprate materials is sup- 
ported by various spectroscopies: spin-resolved photoc- 
mission confirms the singlet character of the first ion- 
ization state, [5] and other probes show that doped holes 
in undoped or lightly-doped cuprates reside primarily on 
oxygens. [6-9] The doping evolution of oxygen content has 
been studied extensively by core-level x-ray absorption 
spectroscopy (XAS). [10-12] In the insulating parent com- 
pounds, early oxygen K-edge XAS experiments (oxygen 
core-electron Is — > 2p transition) show that below the 
main absorption edge a single "pre-peak" exists near 530 
eV,[10] attributed to excitations into the upper Hubbard 
band (UHB). This nonzero projection of the UHB onto 
the oxygen XAS spectra mainly results from a strong hy- 
bridization between copper and oxygen, and the mixing 
of d 9 and d 10 L character in the ground state. Upon hole- 



doping, a lower-energy (d 9 L) peak emerges near 528.5 
eV and grows linearly with hole concentration, while the 
intensity of the higher-energy UHB decreases. Polarized 
x-ray studies indicate that the lower-energy peak in the 
underdoped regime shows almost no oxygen 2p z content, 
but is dominated by planar oxygen 2p XyV contributions, 
[11, 12] lending further credence to the ZRS picture. 

However, recent measurements report that while the 
low-energy ZRS weight increases linearly with doping 
over a wide doping range, it deviates from a linear 
trend and exhibits a weak doping dependence beyond 
~ 20%. [13, 14] It was suggested that this "anomalous" 
behavior implies the inapplicability of the single-band 
Hubbard model and a breakdown of the ZRS picture. [14] 
But cluster dynamical mean-field theory (DMFT) cal- 
culations of the single-band Hubbard model argue that 
the experimental finding could be explained, by showing 
a change in slope with doping of the integrated unoc- 
cupied density of states (DOS) beyond a certain dop- 
ing level. [15] On the other hand, calculations for the 
more sophisticated three-orbital (copper-oxygen) Hub- 
bard model with single-site DMFT find that the low- 
energy weight in the oxygen-projected DOS simply in- 
creases linearly with doping; the result contrasts with 
experiment and challenges even the validity of the three- 
orbital model. [16] Currently, the debate continues and 
questions remain. [14-19] Can oxygen K-edge XAS be 
fully described by the Hubbard models? Are these 
model Hamiltonians and the ZRS picture suitable for the 
cuprates in the overdoped region? 
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FIG. 1. (Color online) (a) Cartoon picture of a Zhang-Rice 
singlet formed by a copper 3d x 2_ y 2 hole hybridized with a 
hole on its neighboring oxygen 2p XtV orbitals. (b) The CusOie 
cluster with periodic boundary conditions considered in the 
exact diagonalization calculations, (c) Schematic of the oxy- 
gen K-edge x-ray absorption process, highlighting additional 
spectral weight reshuffling due to core-hole interaction. 

In this work, we investigate these questions by track- 
ing the doping evolution of photoemission and oxygen 
K-edge XAS spectra in the three-orbital Hubbard model. 
Using large-scale exact diagonalization, we treat many- 
body correlations exactly and include the effect of core- 
hole interaction explicitly, which is essential in the proper 
evaluation of the XAS cross sections. Our calculations 
reproduce the experimental data in a quantitative way 
and explain the crossover or apparent saturation in ZRS 
intensity in the ovcrdopcd regime, highlighting the dy- 
namical spectral weight transfer from the UHB to the 
ZRS band expected for strongly correlated systems. The 
slope of the ZRS weight versus doping further decreases 
due to the presence of core-hole interaction. These re- 
sults confirm the correlated, multi-orbital nature of the 
low-energy quasi-particles in the cuprates, but also in- 
dicate that the ZRS picture may indeed survive to high 
doping levels, albeit with an evolving orbital character. 

Methods- We consider the three-orbital Hubbard model 
on a square lattice of planar Cu02 plaquettes containing 
copper 3d x 2_ y 2 and oxygen 2p x , y orbitals: [20] 

H = - tpdidiaPj* + h.C.) - ]T trf{ P ] aPj , a + h.C.) 

(j,3')a 

ia ja i j 

Here d\ a (p^ a ) creates a hole with spin a at a copper site 
i (oxygen site j), and nf a (n^) is the copper (oxygen) 
hole number operator. The first two terms of the Hamil- 
tonian represent the nearest-neighbor copper-oxygen and 
oxygen-oxygen hybridization. The hopping integrals {t pd 

and t J p J p ) can change sign depending on the phases of the 
overlapping wavefunctions [Fig. 1(a)]. [20, 21] The third 



and fourth terms are the copper and oxygen site energies, 
and the last two terms are the on-site Hubbard interac- 
tions on copper and oxygen, respectively. We use the 
following parameters (in units of eV): e p — = 3.23; 
U d = 8.5, U p = 4.1; \t pd \ = 1-13, \t pp \ = 0.49. These pa- 
rameters have been employed to study the cuprate ma- 
terial La2CuC>4 and produce various spectral features in 
good agreement with experiment. [22, 23] 

Using exact diagonalization for a CusOi6 cluster with 
periodic boundary conditions [Fig. 1(b)], we obtain the 
ground state and calculate the photoemission and XAS 
cross sections for four doping levels: 0% (undoped), 
12.5% (underdoped), 25% (overdoped), and 37.5% (heav- 
ily overdoped). To calculate the spectral functions, we 
exploit translational symmetries of the problem and di- 
agonalize the Hamiltonian matrices in momentum space. 
For oxygen if -edge XAS, the presence of local core- hole 
interaction breaks translational symmetry, and we per- 
form the calculations in real space with a larger matrix 
size. The largest Hilbert space considered in this work 
contains ~ 5.7 x 10 9 basis states. [24] 

Results and Discussion- Oxygen K-edge XAS pro- 
motes Is electrons into empty 2p states [Fig. 1(c)], and 
the resulting spectrum is related to the unoccupied DOS 
measured in inverse photoemission spectroscopy (IPES). 
To capture the doping evolution, we first discuss the 
orbitally- resolved spectral functions [Fig. 2], from which 
we obtain the oxygen-projected DOS. The ground state 
energy of the valence electrons at the corresponding fill- 
ing is defined as zero energy. At 0% doping, the spec- 
tra show an indirect gap of size ^1.5 eV, related to 
charge-transfer excitation of moving a hole from cop- 
per to oxygen. [26] Without doping, the first electron re- 
moval state appears at (ir/2,ir/2), while the first elec- 
tron addition state occurs at (it, 0). This indirect gap 
and its size are consistent with angle-resolved photoemis- 
sion (ARPES) experiments. [27] In an indirect gap sys- 
tem, optical conductivity measurements would obtain a 
gap larger than the actual band gap, as no net momen- 
tum is transferred in the process. We find an optical gap 
~ 1.7 eV, also in agreement with the observed gaps (1.5- 
2.0 eV) in insulating cuprate parent compounds. [28-31] 

Upon hole-doping, low-energy peaks above the Fermi 
level emerge especially at momentum points along the 
magnetic Brillouin zone boundary. This emergent band 
[near zero energy in Figs. 2(b) and 2(c)] consisting of 
both copper and oxygen are associated with the ZRS; its 
orbital composition, spectral intensity, and energy posi- 
tion show strong momentum-dependence. In particular, 
the ZRS band contains vanishing oxygen weight at the T 
point, due to a cancellation in the phase factor of the ZRS 
wavefunction, [32-34] well-defined only away from the T 
point. This is not captured in the single-band model and 
stresses the need for considering the multi-orbital nature 
in describing the ZRS doping evolution. 

On the occupied side [solid lines in Fig. 2], the spec- 
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FIG. 2. (Color online) Spectral functions calculated using 
exact diagonalization on a CusOi6 cluster at (a) half-filling, 
(b) 12.5%, and (c) 25% hole dopings. The ground state energy 
of the valence electrons at the corresponding filling is defined 
as zero energy. The solid and the dashed lines represent the 
electron occupied and unoccupied states, respectively. 

tral features at (tt, 0) and (tt/2, 7r/2) show distinct doping 
dependences: While the weight of the highest occupied 
state at (tt/2, tt/2) gradually decreases upon hole-doping, 
the intensity increases for the highest occupied state at 
(tt, 0). Moreover, the first electron removal state changes 
from momentum (7r/2,7r/2) in the underdoped region to 
(tt, 0) on the overdoped side, reminiscent of the observa- 
tions from ARPES.[35, 36] This doping dependence can- 
not be accounted for by a simple chemical-potential shift 
in a rigid-band model, but is a manifestation of spectral 
weight redistribution due to correlation effects. [37-41] 

The occupied bands between roughly —3 to —5 eV are 
related to oxygen bands dispersing with t pp [42] and the 
Zhang- Rice triplet. [43] In real cuprate materials, other 
orbitals (such as copper 3d,3 Z 2_ r 2 and non-bonding oxy- 
gens) also contribute to the spectral weight in this en- 
ergy range. [44] Below -6 eV [not shown], we find a highly 
incoherent band containing mostly copper character, as- 
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FIG. 3. (Color online) Calculations of oxygen K-edge XAS 
(solid red lines) and IPES measurements of unoccupied 
oxygen-projected DOS (black dotted lines) at various dop- 
ings. The yellow shaded area represents the emergent low- 
energy ZRS band. The XAS and IPES spectra are lined up 
in energy by their peak maximum positions. 

sociated with the lower Hubbard band. These bands at 
deeper binding-energies show a relatively weak doping 
dependence in the doping range considered in this study. 

We next discuss in Fig. 3 the oxygen if -edge XAS cal- 
culations, which include an additional oxygen Is — 2p 
core-hole interaction (%2 i(T(T i U c n]^n^,) with U c = 6 
eV. [45-47] The spectral peaks at 0% doping are closely 
related to the UHB [near 1-2 eV in Fig. 2], which upon 
doping becomes more incoherent and shifts its weight 
to the lower-energy ZRS band. The energy separation 
between these two bands increases systematically with 
doping: the UHB disperses toward higher energy and 
the ZRS band moves in the opposite direction, in agree- 
ment with experiment. [10] Figure 3 also shows the IPES 
calculations of the oxygen unoccupied DOS. The -fT-edge 
XAS lineshape resembles closely that of IPES, because 
the only core-hole interaction is a monopole term form- 
ing a simple charge density attraction due to the isotropic 
Is orbital. [48] However, the two spectra still vary quan- 
titatively [Fig. 4]. [49] 

In Fig. 4 we compare our theory with oxygen K-edge 
XAS experiments on La 2 _ x Sr x Cu0 4 (LSCO).[14] LSCO 
has a relatively simple crystal structure and is suitable 
for systematic studies of electronic properties over a wide 
doping range. As shown in Fig. 4, the experimental ZRS 
weight increases roughly linearly with doping in the un- 
derdoped region. Near optimal doping, the rate of in- 
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FIG. 4. (Color online) Doping dependence of the ZRS spectral 
weight. The La 2 - a; Sr ;c Cu04 (LSCO) XAS data are from Refs. 
10 (Chen), 12 (Pellegrin), and 14 (Peets). The XAS theory 
(red circle) is obtained by integrating the yellow shaded area 
in Fig. 3. The black circle represents the ZRS weight from the 
IPES calculation of the unoccupied DOS. Both theory curves 
are rescaled by normalizing their respective 12.5% ZRS weight 
to the straight line suggested by experiment, [14] and the ZRS 
intensity at 0% doping is defined as zero . 

crease changes and exhibits a weaker doping dependence 
(smaller slope). It was suggested that the ZRS intensity 
may even saturate on the overdoped side (> 20%). [13, 14] 
For comparison, in Fig. 4 we also plot the ZRS weight 
of the XAS and IPES calculations. We define the ZRS 
intensity as zero at 0% doping and rescale both theory 
curves by normalizing their respective 12.5% ZRS weight 
to the straight line suggested by experiment. [14] Neither 
calculation shows a saturation, but only a linear increase 
at low dopings and a change of slope at a higher doping 
level. This crossover in slope is related to the spectral 
weight transfer from the UHB to the ZRS band, which 
can be illustrated by considering the single-band Hub- 
bard model in the atomic limit. [37-41] On a half-filled 
iV-site lattice, the electron removal spectrum (associated 
with singly occupied states) and addition spectrum (as- 
sociated with the unoccupied UHB) both have a spec- 
tral weight equal to N. Upon doping x holes, there are 
N — x singly occupied states, and the UHB weight also 
becomes N — x. Due to the conserved sum of the occu- 
pied and unoccupied states (2N), the low-energy spectral 
weight which emerges upon doping x holes is thereby in- 
creased by 2x. In contrast, for an uncorrelated system 
the emergent spectral weight equals x when doped with 
x holes. The occurrence of spectral weight transfer sig- 



nals the presence of strong correlations, and the effect 
weakens when the system becomes less correlated. 

Although we do not find an apparent saturation of the 
ZRS weight with doping, our XAS calculation seems to 
capture well the experimental findings. In particular, 
the XAS theory curve matches the experimental data 
within error bars near optimal doping and in the over- 
doped regime; our theory also predicts the doping de- 
pendence of the heavily overdoped (x ~ 0.35) LSCO 
sample. Compared with IPES (U c = 0), the core- hole 
attraction (U c — 6 eV) effectively overcomes the energy 
cost of double occupancy and shifts roughly an additional 
~ 10% of the UHB weight to lower energies. This ef- 
fect is more prominent at low doping where a substantial 
UHB weight still remains, leading to a more pronounced 
change of slope in ZRS intensity in XAS. The IPES cal- 
culations underestimate the low-energy spectral weight at 
low dopings and overestimate the slope of the ZRS inten- 
sity with doping in the overdoped region. We note that 
Liebsch[17] and also Peets et al. [19] have pointed out the 
importance of choosing a proper energy integration win- 
dow for the ZRS band when comparing theory and ex- 
periment. Here we use a ~2 eV energy window whose 
upper limit changes with doping and resides in the min- 
imum between the centroids of the UHB and the ZRS 
band (see Fig. 3), consistent with experiment. [14, 19] 

As mentioned previously, the doping evolution of the 
low-energy ZRS spectral weight has also been studied by 
various theoretical techniques. [15-17, 39, 50, 51] In par- 
ticular, single-site DMFT using the three-orbital model 
shows a linear increase with doping of the low-energy 
weight. [16] Cluster DMFT of the single-band model 
demonstrates a change in slope of the ZRS weight be- 
yond a certain doping level, [15, 17] which resembles the 
behavior of our calculations for the unoccupied DOS in 
the three-orbital model. However, the slope on the over- 
doped side is expected to be even larger in the single-band 
case, [40] as there the UHB completely disappears close 
to 100% hole doping, while in the three-orbital model the 
d 8 double occupancy, albeit small, stays nonzero. 

To summarize, our exact diagonalization calculation 
of the three-orbital Hubbard model agrees with the XAS 
experiments at a quantitative level, which relies crucially 
on the combined effect of correlation-induced spectral 
weight transfer and an additional spectral weight redis- 
tribution due to core-hole interaction. The change of 
slope in ZRS intensity as a function of doping manifests 
clearly the presence of strong electronic correlation. We 
also find a strong doping and momentum dependence of 
the orbital composition in the ZRS band, which cannot 
be fully described by the single-band Hubbard model. 
Out results potentially confirm the continuation of the 
ZRS picture even in the overdoped region. Future ex- 
periments on overdoped samples would be helpful to il- 
luminate the nature of the emergent low-energy quasi- 
particles in cuprate superconductors. 



■5 



Acknowledgments- This work is supported by the U.S. 
DOE under Contract No. DE-AC02-76SF00515. C.C.C. 
is supported by the Aneesur Rahman Fellowship at ANL. 
Y.F.K. is supported by the U.S. DOD through the ND- 
SEG Fellowship. C.J.J, is supported by the Stanford 
Graduate Fellowships. R.T. is supported by an SITP Fel- 
lowship at Stanford University and SPP 1458/1. A.P.K. 
acknowledges support from the DFG through TRR 80. 
The simulations were performed on the Hopper peta-flop 
Cray XE6 system at NERSC, supported by the U.S. DOE 
under Contract No. DE-AC02-05CH11231. 



[1] 
[2] 

[3] 

[4] 

[5] 
[6] 
[7] 
[8] 

[9: 

[io; 
[ii 

[12 
[13 
[14 
[15 

[16 

[17 
[18 
[19 
[20 
[21 
[22 
[23 
[24 



J. Hubbard, Proc. R. Soc. Lond. A 276, 238 (1963). 
K. A. Chao, J. Spalek, and A. M. Oles, Phys. Rev. B 18, 
3453 (1978). 

F. C. Zhang and T. M. Rice, Phys. Rev. B 37, 3759(R) 
(1988). 

H. Eskes and G. A. Sawatzky, Phys. Rev. Lett. 61, 1415 
(1988). 

N. B. Brookes et al, Phys. Rev. Lett. 87, 237003 (2001). 

Z.-X. Shen et al, Phys. Rev. B 36, 8414 (1987). 

J. M. Tranquada et al, Phys. Rev. B 35, 7187 (1987). 

N. Niicker et al, Phys. Rev. B 37, 5158 (1988). 

Y. Sakurai et al, Science 332, 698 (2011). 

C. T. Chen et al, Phys. Rev. Lett. 66, 104 (1991). 

C. T. Chen et al, Phys. Rev. Lett. 68, 2543 (1992). 
E. Pellegrin et al, Phys. Rev. B 47 3354 (1993). 

M. Schneider et al, Phys. Rev. B 72 014504 (2005). 

D. C. Peets et al, Phys. Rev. Lett. 103 087402 (2009). 
P. Phillips and M. Jarrell, Phys. Rev. Lett. 105 199701 
(2010). 

X. Wang, L. de' Medici, and A. J. Millis, Phys. Rev. B 
81, 094522 (2010). 

A. Liebsch, Phys. Rev. B 81, 235133 (2010). 

D. C. Peets et al, Phys. Rev. Lett. 105 199702 (2010). 

D. C. Peets et al, arXiv: 1004. 1146. 

V. J. Emery, Phys. Rev. Lett. 58, 2794 (1987). 

E. Arrigoni et al., New J. Phys. 11, 055066 (2009). 
Y. Ohta et al, Phys. Rev. B 43, 2968 (1991). 

C.-C. Chen et al, Phys. Rev. Lett. 105, 177401 (2010). 
To obtain the ground state, we use both the Lanc- 
zos tridiagonalization technique and the Parallel ARnold 
PACKage based on implicitly restarted methods. [25] We 
adopt a block-checkerboard matrix decomposition with a 
parallel computing scheme implemented on a distributed 
memory system. By partitioning both the rows and 



columns of a matrix, each computing thread no longer 
stores a full vector, and the matrix diagonalization is not 
limited by the total memory per computing node. 

[25] R. B. Lehoucq, D. C. Sorensen, and C. Yang, ARPACK 
Users' Guide: Solution of Large-Scale Eigenvalue Prob- 
lems with Implicitly Restarted Arnoldi Methods (Philadel- 
phia: SIAM, 1998). 

[26] J. Zaanen, G. A. Sawatzky, and J. W. Allen, Phys. Rev. 
Lett. 55, 418 (1985). 

[27] N. P. Armitage et al, Phys. Rev. Lett. 88, 257001 (2002). 

[28] Y. Tokura et al, Phys. Rev. B 41, 11657 (1990). 

[29] S. L. Cooper et al, Phys. Rev. B 42, 10785 (1990). 

[30] S. Uchida et al, Phys. Rev. B 43, 7942 (1991). 

[31] J. P. Falck et al, Phys. Rev. Lett. 69, 1109 (1992). 

[32] P. Unger and P. Fulde, Phys. Rev. B 48, 16607 (1993). 

[33] J. J. M. Pothuizen et al, Phys. Rev. Lett. 78, 717 (1997). 

[34] Q. Yin et al, Phys. Rev. Lett. 100, 066406 (2008). 

[35] C. Kim et al, Phys. Rev. Lett. 80, 4245 (1998). 

[36] A. Ino et al, Phys. Rev. B 62, 4137 (2000). 

[37] E. Dagotto et al, Phys. Rev. Lett. 67, 1918 (1991). 

[38] H. Eskes,, M. B. J. Meinders, and G. A. Sawatzky, Phys. 
Rev. Lett. 67, 1035 (1991). 

[39] M. S. Hybertsen et al, Phys. Rev. B 45, 10032 (1992). 

[40] M. B. J. Meinders, H. Eskes, and G. A. Sawatzky, Phys. 
Rev. B 48, 3916 (1993). 

[41] M. A. van Veenendaal et al, Phys. Rev. B 49, 1407 
(1994). 

[42] Y. Mizuno et al, Phys. Rev. B 57, 5326 (1998). 
[43] T. Learmonth et al, Europhys. Lett. 79, 47012 (2007). 
[44] C. Weber, K. Haule, and G. Kotliar, Phys. Rev. B 82, 
125107 (2010). 

[45] In copper if-edge XAS, the Is core-hole interaction is 
found to be slightly larger than the on-site Hubbard in- 
teraction based on the resonant energies of the well- and 
poorly-screened states. For oxygen K-e&ge XAS, the oxy- 
gen Hubbard repulsion U p is estimated to be 4-5 eV, and 
a core-hole strength U c ~5-6 eV is commonly used. 

[46] H. Eskes and G. A. Sawatzky, Phys. Rev. B 43, 119 
(1991). 

[47] K. Okada and A. Kotani, J. Phys. Soc. Jpn. 76, 123706 
(2007). 

[48] In transition- metal L-edge XAS (2p — > 3d transition), 
core-hole interactions involve high-order multipole terms 
in the expansion of the orbital angular wavefunctions. 
This leads to XAS spectra with a complicated multiplct 
structure whose lineshapes and spectral weights differ 
substantially from the original unoccupied DOS. 

[49] C. J. Jia et al, arXiv: 1109.3446. 

[50] K. Okada and A. Kotani, J. Phys. Soc. Jpn. 75, 123703 
(2006). 

[51] T. Ahmed et al, Phys. Rev. B 83, 115117 (2011). 



